Fast Computation of First-Order Feature-Bispectrum Corrections 
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Features in the inflaton potential that are traversed in much less than an e-fold of the expansion can 
produce observably large non-Gaussianity. In these models first order corrections to the curvature 
mode function evolution induce effects at second order in the slow roll parameters that are generically 
greater than ~ 10% and can reach order unity for order unity power spectrum features. From a 
complete first order expression in generalized slow-roll, we devise a computationally efficient method 
that is as simple to evaluate as the leading order one and implements consistency relations in a 
controlled fashion. This expression matches direct numerical computation for step potential models 
of the dominant bispectrum configurations to better than 1% when features are small and 10% when 
features are order unity. 



I. INTRODUCTION 

Features in the inflaton potential can give rise to large 
non-Gaussianity [HE]. In order to also satisfy cosmic 
microwave background (CMB) constraints on the power 
spectrum, observability in the bispectrum of an individ- 
ual feature requires that it be traversed when the hori- 
zon was of order the current horizon in a small fraction 
of an e-fold [3]. As this represents a strong violation of 
the slow-roll limit, corrections to the leading-order bis- 
pectrum expression can approach order unity. Curiously 
such models are in fact favored by the WMAP 7-year 
CMB power spectrum data and so further tests require 
more accurate techniques [3j. The improvement of the 
fit due to these models is not due to the glitch in the 
CMB data at ^ = 20 — 40 "^-"TT , but comes predominantly 
from the region around first peak. Similar improvements 
have also been noted arising from transplanckian mod- 
ifications to the power spectrum f8], axion monodromy 
inflation [9 and non-trivial initial vacuum choices |10j . 
While the polarization power spectrum is likely to pro- 
vide the next test of such models HI] , the bispectrum 
of these models may also provide an important consis- 
tency check. 

Features in the inflationary potential have a long his- 
tory and have been considered by many authors going 
back to Starobinsky [T^] who first calculated the power 
spectrum due to a potential with a discontinuous first 
derivative. The bispectrum due to such a feature was 
calculated by jl3], see also [M]. The fiuctuations in 
the power spectrum due to a step feature in the po- 
tential were considered by [THl [H] and their resulting 
non-Gaussian signatures in the bispectrum were calcu- 
lated numerically by [H [2] and an analytic approxima- 
tion was found by [31 [T7] . Steps in the warp factor and 
potential in the context of brane infiation were consid- 
ered by |18j . Features in the power spectrum and large 
non-Gaussianities arising from resonant effects in the in- 
flationary potential were first proposed by j2] and then 
discovered in axion monodromy inflation [^ [THl [5D] . 

An exact calculation of the bispectrum for such models 
requires a computationally intensive direct integration of 



the curvature fluctuations for each Fourier space config- 
uration [3 12] and is impractical for analysis purposes. 
In [17], a fast approximate method was developed based 
on the generalized slow roll (GSR) approach. Here the 
curvature fluctuation mode functions are iteratively cor- 
rected due to the presence of the feature. The leading- 
order expression was shown to be accurate to typically 
ten percent for small amplitude features and up to order 
unity for large amplitude features. In the small ampli- 
tude feature limit, this error is associated with slow-roll 
corrections that change the phase of the modefunctions 
[3]. Even without a feature, these corrections typically 
generate 10% changes to the bispectrum [3T]. A full cal- 
culation of all first-order correction terms is computation- 
ally cumbersome |17j and again becomes impractical. 

In this paper we show that the dominant first-order 
correction can be simply computed. In pi] we review the 
GSR approach and give all next to leading-order terms in 
the power spectrum and bispectrum. In §1111 we isolate 
the dominant terms, test them against exact calculations, 
and establish their compatibility with power spectrum 
corrections in the squeezed limit. We discuss these results 
infJIVl 



II. GSR APPROXIMATION 

In this section we review the GSR formalism [THJ [H] 
for computing the power spectrum ^23' and bispectrum 
|17) for inflationary models with relatively large ampli- 
tude features including next to leading-order corrections. 
Such corrections are flrst-order in the GSR iteration and 
second order in the deviations from slow-roll. 



A. Power Spectrum 

Beyond the slow-roll approximation, the curvature 
power spectrum can be computed exactly in linear theory 
as 
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where x = krj and rj is the conformal time to the end of 
inflation. The modefunction y satisfies the Mukhanov- 
Sasaki equation, 
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which is not necessarily small or constant. Throughout 
we set the reduced Planck mass Mpi = (SttG)"^/^ = 1 
as well as ?i = c = 1. Note that /"// in the source 
function g involves first derivatives of the second slow 
roll parameter 
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Up to this point, no assumptions have been made be- 
yond the validity of linear theory but the modefunction 
y remains an implicit functional of en- Briefly, the GSR 
approach to solving the modefunction equation ([2| is to 
consider the RHS as an external source with an iterative 
correction to y [TB]. To lowest order, we replace y ^ yo 
where 
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is the solution to the equation with 5 — > 0. The first-order 
mode function y can then be obtained through the Green 
function technique. The result is an approximation that 
still requires ch to be small in an absolute sense but 
allows it to evolve rapidly with large fractional changes 
to its value. Hence \g\ can be large in an absolute sense. 
To leading order, the curvature power spectrum is 
given by 
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where kr]^ <^ 1 and 



^j^, , 3sin(2ii) 3cos(2m) 3sin(2u) 

^m^^TTj, -3^ ^- — • (9) 



2u3 



2m 



The source function for the power spectrum is given by 



G = -21n/+-(ln/)', 



(10) 



and thus 

G' = -2(ln/)' + ^(In/)" = ^g - g[(ln/)f . (11) 

Note that the leading order expression for the power spec- 
trum is already first-order in the deviations from slow 
roll. 

The addition of the term quadratic in (In /)' to g in 
Eq. ( 11 1 makes G' a total derivative and guarantees that 



the power spectrum is independent of the arbitrary epoch 
r]^ after horizon crossing, ensuring that the curvature re- 
mains constant thereafter |23) . 

To first order in the mode function iteration and second 
order in the slow-roll parameters [22l [23] 
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When /"// controls the large deviations in G', the dom- 
inant term is /i 23j 
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For a wide range of models, this has been shown to be a 
good approximation when |/i| < l/-\/2 [24]. In this ap- 
proximation, the power spectrum depends only on a sin- 
gle source function G' through two single integrals over 
the functions W and X. We seek a similarly simple but 
accurate approximation for the bispectrum in what fol- 
lows. 



B. Bispectrum 

For models in which large slow-roll corrections arise 
from a sharp feature where 77^ or /"// becomes large, 
the bispectrum can be approximated as |17) 

Bn{h{) = ml inkAv*)TZkAv*)T^kAv*) (17) 
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where the curvature fluctuation is given by 
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the shorthand convention k^ =^ ki,k2,k^, and 77* is an 
epoch well after all modes have left the horizon. 

Using the same iterative GSR approach to evaluate y, 
the bispectrum to first-order in the modefunction correc- 
tion or second order in the slow-roll violation becomes 
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The leading order response to this source comes from 
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The first-order correction terms again involve G' and 



UiAikii]) 



V2 ydlnri 
X '^[yo{kiv)yo{k2V)[yoihv) + yoiksij)]] 
+ cyc., 

UMhT^) =-\ r ^G'{lnf,)W{ksfi) 

27,, V 
X (j^ + A ^{yoikiv)yo{k2v)yo{k3v)] 



UicikiT]) 



+ cyc., 

/i(fc3) , 1 '•°° 
2 



V2 
d 

din 77 
eye, 
3 



^G'{lnfj)X{k3fj) 
n ''I 



3 '^[yaikir])yo{k2V)yoik3V)] 



1 



1 



UMhv)=-j I ^G'(lnry) , . 

X f ^1^ + 3J 3[yo(fci77)yo(fc2??) 

[yolks'?) + yo(fc3'7)]]+ eye, 
UiEihr]) = - 33?[2/o(fci'?)yo(fc2?7)yo(fc3?7)] 



/ 



1 //' 



2ffs V/ 



(22) 



where eye. denotes the 2 additional cyclic permutations 
of the k indices. This first-order expansion has been 
shown to be highly accurate even for large amplitude 
features but has the drawback that it is cumbersome to 
compute since the nested integrals in Uib-d involve con- 
figuration dependent quantities |17j . 



III. FAST BISPECTRUM COMPUTATION 

While first-order corrections to the bispectrum for 
models with features are generally at least of order 10% 
and hence important for accurate computation, most of 
the first-order terms in Eq. ( 22 ) arc irrelevant where the 

On the other hand some 



bispectrum is observably large, 
of these terms are important for maintaining physicality 
in the superhorizon and squeezed limits. In |III A| and 
IIIB we devise and then test a first-order methodology 



that both efficiently corrects the bispectrum where it is 
large and implements physicality constraints in a con- 
trolled fashion. 



A. Methodology 

An observably large bispectrum arises if potential fea- 
tures are sufficiently sharp that they are traversed in 
much less than a Hubble time. In this case, the dominant 
contributions to the bispectrum arise when the modes 
are well inside the horizon krj ^ \. Even if the features 
themselves are of small amplitude, first-order corrections 
from the slow-roll contributions to G' make substantial 
fractional corrections to the bispectrum [3]. In this case 
G" is nearly constant and we can evaluate the expressions 
in Eq. (22 1. We derive our method from these two as- 



sumptions but show that the resulting expressions remain 
a good approximation even when features in the power 
spectrum reach order unity or when one of the sides of 
the triangle is outside the horizon. 

Under these two assumptions, all the nested integrals 
contribute very little to the integrand in Eq. ( 19 ) in that 



most of their impact is around or after horizon crossing 
rather than before. In addition, models with large bispec- 
tra contribute mainly through the /"// terms in gs and 
G' and so Uie is negligible. The remaining terms are the 
ones proportional to Ii{k) in Uia and Uic- In fact the 
yoyayo contributions of these two terms cancel exactly 
leaving only yoyoyo contribution similar to the zeroth or- 
der term (21). As noted in [5], this fact means that for 
triangles where all three kr] ^ 1, the first-order correc- 
tions can be cast in the same form as that of the zeroth 
order term, namely in terms of single integrals involving 
only the perimeter of the triangle K = ki + k2 + k^. 

As kr] decreases, the Ji and integral contributions to 
Uic cancel, leaving the Uia term as the dominant cor- 
rection. However as ki] — >■ 0, the other first-order correc- 
tions ensure that the bispectra remain constant in accord 
with conservation of the comoving curvature TZk- This 
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FIG. 1. Equilateral bispectrum for a small amplitude step c = 10^'' with d = 0.0001 (left) and 0.0003 (right). Upper panel 
shows the GSRi first order approximation versus the full numerical computation. Lower panel shows the percentage difference 
between the two as well as between the zeroth order approximation GSRo and the numerical one. 



condition can be alternately enforced by replacing the 
bispectrum source with a total derivative |17) 



and 



.9B^G'B(ln77) 



Ah -VhV 



(23) 



V / 

just as in the case of the power spectrum. Combining 
these considerations, our fast first-order bispectrum ex- 
pression becomes 
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WBoix) = xsina::, 
WBiix) — cos a;, 
WB2ix) — sinx/a;, 



Wb3{x) = — xcosx, 
Wsiix) = sinx, 
Wb5{x) = — cosx/a:. 



The lack of a. K — 2kk term in Ibs comes from the a; ^ 1 
cancellation of the Uia and Uic terms. 

We shall call the n = — 2 terms the zeroth order 
approximation and the full set n = — 5 the first or- 
der approximation, denoting these GSRq and GSRi re- 
spectively. Note that the boundary term at Krj^ <C 1 
is formally only defined correctly for the n = 1,2 terms 
where Wbu — > 1 but we retain the others for notational 
compactness as their contributions vanish in this limit. 
While the Ib5 term diverges as Krj^ — >• we shall see 



that the differencing construction in Eq. ( 24 ) guarantees 



that this divergence or the corresponding dependence on 
the arbitrary scale ry* has no observable consequence. 

The computational cost of GSRi is only double that of 
GSRo involving 6 rather than 3 one dimensional integrals 
in addition to the power spectrum correction Ii{k) from 
Eq. ( [T3| . That the corrections are all proportional to 
Ii{k) is an important feature in our construction and 
we shall see enforces compatibility between the power 
spectrum and bispectrum approximations. 

Although this construction is motivated by the slow- 
roll corrections to the bispectrum, we shall see next that 
it also works quite well for case where the correction is 
dominated by the feature itself. 




FIG. 2. Equilateral bispectrum for a large amplitude step c = 5.75 x 10^* with d = 0.0003 (left) and 0.0001 (right). Upper 
panels shows the GSRi first order approximation versus the full numerical computation. Lower panels shows the percentage 
difference between the two as well as between the zeroth order approximation GSRo and the numerical one. 



B. Configuration Tests 

As an example of a feature model with a possibly ob- 
servable bispectrum, we consider the potential |15| 



V{<P) = ^m^^ 
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which corresponds to a smooth step ai (j) = (ps of frac- 
tional height c and width d. Given the WMAP7 pref- 
erence for a large scale feature, we set (jjs to repro- 
duce the scale of the maximum likelihood (ML) solution 
ri{4>s) = 8.163 Gpc [3 . For the amplitude, we take the 
case where c — >■ and the maximum likelihood solution 
for the WMAP7 data c = 5.75 x 10"^. For the width d 
we take several cases that would lead to observable bis- 
pectra with this amplitude and position. We choose m 
to be compatible with the WMAP7 slope for the under- 
lying slow- roll potential m = 7.126 x 10~^ such that the 
underlying tilt from the smooth part of the potential is 
fis « 0.963. Following the notation of the existing litera- 
ture [11 [2j , we construct the bispectrum statistic 
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where we take As = 2.39 x 10~^ which is approximately 
the amplitude of the power spectrum in the absence of 
the feature. 

Most of the observable impact of features in the bispec- 
trum comes from equilateral triangles where fci ^ fc2 ~ k^ 
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and so as fcyy — ^ the quantity in brackets vanishes and 
the expression becomes independent of the arbitrary end 
point ry*. 

In Fig. n] (upper), we compare the first order approxi- 
mation (GSRi) to the full numerical calculation for the 
equilateral case and a small amplitude step c = 10^^ with 
d = 0.0001 and 0.0003. Differences in the upper panel 
are quoted as percentages of a smooth envelope 
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in the lower panel to avoid dividing by an oscillatory 
quantity. Here the equilateral damping scale 
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and eo = 0.00925 is the value calculated on the slow- 
roll background in the absence of the step. The form of 
the envelope is derived from the analytic solutions in [3] 
and is accurate at the several percent level for all models 
considered here. 

For c = 10^^, the agreement of GSRi and the exact 
numerical treatment is excellent, with differences below 
the 1% level between the first oscillation and the damp- 
ing scale set by the finite width d (upper panel) . The first 
order correction eliminates the f 0% errors of the GSRq 
approximation that appear mainly as a phase error (lower 
panel). Deep in the damping tail, the first order solution 
develops a small phase difference leading to larger frac- 
tional errors but controlled amplitude errors. This error 
can be attributed to terms in the first order expansion 
at Eqn. ( 22 ) that are not captured by the approximation 
at Eqn. (24). In particular, contributions due to some 



of the nested integrals are damped more slowly than the 
leading order contributions. This means that deep in the 
damping tail they can become a significant fraction of 
the leading order result. However, since this effect only 
becomes important in the region where the bispectrum is 
already small, we conclude that it can be safely ignored. 
This good agreement persists until the feature makes 
order unity changes to the power spectrum. In Fig.[2J we 
show the maximum likelihood amplitude model. Here the 
oscillations take on a modulated form where every third 
extrema is reduced in amplitude reflecting the strong os- 
cillations in the power spectrum. These are well captured 
by the first order approximation between the first oscil- 
lation and the damping tail with residuals around ^ 6% 
correcting the ^ 40% errors of the zeroth order approx- 
imation. Errors in the damping tail grow again mainly 
due to a phase error but remain small until the bispec- 
trum amplitude has damped to an unimportant level. 
Note that in this example the maximum of |/i| w 0.25. 
Like the power spectrum, this quantity monitors the ac- 
curacy of the first order computation. The criteria for 
equilateral bispectrum is slightly more stringent than the 
power spectrum due to the 3 fc-modes that can be cor- 
rected and hence 
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FIG. 3. Squeezed limit bispectrum for a step witli small 
amplitude c — 10^'', d = 0.0003. Shown here is 
12/NL(fcL, A;_L, fcs)/5 for various values of the short side, ks- 
In the upper panel, we show the result of evaluating the first 
order approximation GSRi versus the full numerical compu- 
tation. The lower pa nels here show the slow roll ks rescaling 
GSRis from Eq. (itsl. 



the leading order corrections when all three modes are 
subhorizon scale during feature crossing. 

For squeezed configurations, our first-order approxima- 
tion (24) becomes 
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is the rough criteria for better than 10% percent level 
accuracy. 

An important check of the physicality of the bispec- 
trum is the squeezed limit where ki — ks ^ k2 ~ k'^ = 
k^ and the the bispectrum must satisfy consistency with 
the power spectrum [5S] 
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Note that this is a non-trivial check on our construction 
even for small amplitude features since we retained only 
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the expression becomes independent on the arbitrary end 
point ?7*. 

To check the consistency relation, it is useful to com- 
bine all of the terms to form a single integral over the 
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FIG. 4. Squeezed limit bispectrum for a step with large ampli- 
tude c = 5.75 X lO""* and d = 0.0003. Shown here is the GSRi 
and slow-roll corrected GSRis approximations compared with 
the numerical computation as in Fig. [3] Note that the latter 
corrects all of the ks dependent error. 
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with /i = Ii^k^). Using the approximation of Eq. (16) 
in the I III < 1 limit 
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To establish the consistency relation, we need to relate 
Wsq to Wp and G'g to G' . We can manipulate the latter 
pair via integration by parts 
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For bispectra that are dominated by f"/f, Gb ~ 
— {f'/P)- If we further take the approximation that 
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const., the consistency relation is satisfied 
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Thus the consistency relation holds in our first order bis- 
pectrum approximation as long as / remains nearly con- 
stant. 

Of course / remains constant only at zeroth order in 
slow roll so there are additional correction in a full first 
order calculation. In Fig. p^ (upper), we show that corre- 
spondingly there is a small amplitude mismatch between 
the numerical and the GSRi results which grows loga- 
rithmically with decreasing ks- Recall that the quan- 
tity f^Nhiks, k^^ki,) should become independent of ks as 
ks —>■ to satisfy the consistency relation. 

One can gain further insight into the missing term by 
re-examining the slow-roll limit of the full first-order ex- 
pression in Eq. ( 22 ) beyond the subhorizon approxima- 
tion. Here G" w 1 — n^, where recall fis is the tilt that 
describes the power spectrum in the absence of the slow- 
roll violating feature. The leading contributions are from 
the UiB and Uid terms and account for the evolution of 
/ between horizon crossing 77 « t/ks, when the curvature 
fluctuation for ks froze out, and 77/, the epoch of slow- 
roll violation. They form a multiplicative correction to 
hs and I2B of 



R: 



{' 



^In 
1, 



ksVf > a^max, 



(43) 



with Xniax = e '''^/2 w 2.07 where je is the 
Euler-Mascheroni constant. Note that the ln(fc?7/) 
term is simply the leading order slow-roll expansion of 
ATi{ks)/ATi{r]J ) as one might expect. Furthermore, 
the GSR technique can be shown to imply a slow-roll 
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/2 rather than — 1 (see |16j . 



Eq. 105). In Fig. Is] (lower), we demonstrate the effect of 
this rescaling with Eq. (43), denoted GSRis- Here the 



^5 dependent discrepancy disappears entirely. 

Applying this correction does not impact any triangle 
where all three modes are subhorizon at rjf and hence has 
very little impact on high signal-to-noise modes. Further- 
more given that rjf must be comparable to the horizon 
size for the bispectrum features to be detectable in the 
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FIG. 5. Consistency relation test for a large amplitude step 
c = 5.75 X 10"* and d = 0.0003 (upper) and 0.0001 (lower). 
Here the change in the power spectrum slope due to the fea- 
ture fis — Us and the squeezed bispectrum from the feature 
/nl are both evaluated numerically. For the latter we take 
ksVf ~ 1-38 X 10"^, though the results are independent of 
this choice. 



CMB, observable triangles cannot acquire large logarith- 
mic corrections. We conclude that for practical purposes, 
this correction can be safely ignored. 

At higher values of the step height, c, our scaling still 
removes the fcg dependent errors (see Fig. [4]). This is be- 
cause these errors correspond to the slow evolution of / 
between rj — 1/ks and rjf rather than the feature itself. 
On the other hand, a new scaling offset develops that can 
be attributed to the change the feature makes on / at 77/ . 
In principle these could be corrected by integrating the 
feature contributions to G" in the first order Uib and Uijj 
terms. However this offset is of the same order and na- 
ture as those found for the equilateral cases. Given that 
squeezed triangles do not dominate the signal-to-noise, 
we again conclude that further correction is unnecessary. 

In Fig. [5] as a test of the consistency relation itself 
for a large amplitude step, we compare the numerical 
power spectrum slope to the squeezed bispectrum. The 
operator which dominates the bispectrum in the case of 
a feature is a subleading contribution to the usual slow- 
roll consistency relation. Thus, we subtract off the slow- 
roll contribution to the slope of the numerical power 
spectrum, and instead compare the squeezed limit bis- 
pectrum with the deviation away from the slow roll re- 
sult for the slope of the power spectrum, fis — ns(kL), 
where n^ « 0.963 is the slope of the power spectrum 
in the absence of the feature. In practice this correc- 



FIG. 6. Flat limit bispectrum for a large amplitude step 
c — 5.75 X 10""' with d = 0.0003. Flat triangle agreement 
between the GSRi approximation and numerical evaluation 
is comparable to equilateral triangles. 



tion is negligible if kj^rjf ^ 1 or more generally when 
I/atlI ^ C(1 — ns). It is only that the agreement is so 
good that we take it into account here. Furthermore, the 
excellent agreement checks the accuracy of the numerical 
bispectrum calculation. 

Finally it is interesting to note that flat configurations 
where ki — 2fc2 — 2/03 — kp are somewhat special in 
Eq. (24), since the argument to the 75 integrals vanishes 



leaving the expression ill-defined. It is in fact well-defined 
if we take instead the limit of a flat triangle k^ — K — 
2kF -^0 



BnikF,kF/2,kF/2) 



1,6 



AAnikF)^l{kF/2) 



{ - /Bo(2fcF) - 7lBli2kF) + l2lB2{2kF) 

, h{kF 



V2 



-Gj^iBbih 

kp 
2h{kp/2) 



x/2 



lB:i{2kF) - 7lBii2kF) + l2lB^{2kF) 

(44) 

iBsi^kp) - llBii'ikF) + 2lBi{kF) 



-12/B5(2fcp)-6/B5(fcF)]} 



which again becomes independent of 77* as the triangle 
flattens. In Fig.|6J we compare our first order approxima- 
tion to the numerical results. For scales between the first 
oscillation and the damping tail, the approximation cap- 



tures the flat and equilateral bispectrum behavior com- 
parably well. 



IV. DISCUSSION 

We have devised an efficient method to compute first- 
order corrections to the bispectrum of inflationary models 
with features. First order corrections to the bispectrum 
are generically at least ^-^10% and can reach order unity 
for order unity features in the power spectrum. 

Based on the GSR approach, we have shown that cor- 
rections in the high k limit, where the bispectrum ampli- 
tude is set well within the horizon, take on a simple form 
involving only single integrals over the slow-roll param- 
eters. This limit is observationally important since it is 
here that the bispectrum can become large if the inflaton 
crosses a feature in the potential in much less than an 
e-fold of the expansion. We have constructed these cor- 
rections such that they implement consistency relations 
at low k in a controlled fashion. 



Comparison with direct numerical computation of the 
bispectrum shows that the approximation works ex- 
tremely well for the full range where the bispectrum is 
large. For cases where the zeroth order expressions devi- 
ate by 50%, the first order approximation deviates by less 
than 10%. These techniques should be useful for analyz- 
ing any model where inflaton potential features provide 
observably large non-Gaussianity. 
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